Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 4 May 2012 (MN KTeX style file v2.2) 



The intergalactic medium thermal history at redshift z = 1.7-3.2 
from the Ly-a forest: a comparison of measurements using wavelets 
and the flux distribution 



o- 

■ A. Garzillii*, J.S. Bolton^, T.-S. KW, S. Leach\ M. Viel^'^ 

^ SISSA, Astrophysics Sector, via Bonoinea 265, 34136 Trieste, Italy 

_ ^ School of Physics, University of Melbourne, Parkville, VIC 3010, Australia 

^5 • ^ Department of Astronomy, University of Wisconsin-Madison, 473 N. Charter St., Madison, WI53706. USA 

■ ^ INAF-Osservatorio Astronomico di Trieste, via G. B. Tiepolo II, I-34I3I Trieste, Italy 
Cf^ ' ^ INFN/National Institute for Nuclear Physics, via Valerio 2, 1-34I27 Trieste, Italy 



Accepted — . Received — ; in original form 4 May 2012 



ABSTRACT 

We investigate the thermal history of the intergalactic medium (IGM) in the redshift inter- 
val z = 1.7-3.2 by studying the small-scale fluctuations in the Lyman-a forest transmitted 
flux. We apply a wavelet filtering technique to eighteen high resolution quasar spectra ob- 
tained with the Ultraviolet and Visual Echelle Spectrograph (UVES), and compare these data 
to synthetic spectra drawn from a suite of hydrodynamical simulations in which the IGM 
thermal state and cosmological parameters are varied. From the wavelet analysis we obtain 
estimates of the IGM thermal state that are in good agreement with other recent, independent 
wavelet-based measurements. We also perform a reanalysis of the same data set using the 
Lyman-a forest flux probability distribution function (PDF), which has previously been used 
to measure the IGM temperature-density relation. This provides an important consistency test 
for measurements of the IGM thermal state, as it enables a direct comparison of the con- 
straints obtained using these two different methodologies. We find the constraints obtained 
from wavelets and the flux PDF are formally consistent with each other, although in agree- 
ment with previous studies, the flux PDF constraints favour an isothermal or inverted IGM 
temperature-density relation. We also perform a joint analysis by combining our wavelet and 
flux PDF measurements, constraining the IGM thermal state at z = 2.1 to have a temperature 
at mean density of Tq/[10'^ K] = 17.3 ±1.9 and a power-law temperature-density relation ex- 
ponent 7 = 1.1±0.1(1(t). Our results are consistent with previous observations that indicate 
there may be additional sources of heating in the IGM at z < 4. 
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1 INTRODUCTION 

The intergalactic medium (IGM) is the largest reservoir of baryonic 
matter in the early Universe, and so gaining an understanding of its 
physical state and chemical composition is an important goal of 
modern cosmology. In the current picture for the evolution of the 
baryons, there are two reionisation events which turned the neu- 
tral gas in the IGM into an ionised medium. The first reionisation 
event is thought to be caused by hydrogen (and neutral helium) 
ionising radiation produced by early galaxies. The precise redshift 
of this reionisation event is not well constrain ed but it is though t 
to initiate at a redshift no later than z — II jUarson et al.ll201 ih 
and end by 2 ~ 6, which is when the Universe becomes trans- 
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parent to redshifted Lyman-a photons from quasars jSecker et al.l 
2001; Fan et al. 2006). The second reionisation event is expected 
to instead be driven by quasars at lower redshifts, which pro- 
duce a hard ionising spectrum that can reionise singly ionised 
helium hy z 3 jMadau et al.1 1 19991 : iFurlanetto & oi] I2OO8I: 
iMcOuiim et alj|2009f) . Photo-heating during both of these reioni- 
sation events leaves a 'footprint' on the thermal state of the IGM; 
determining the redshift evolution of the IGM temperature can 
therefore he lp pin do wn the details of these reionisation eras (e.g. 
iTheuns et al . 2002; H ui & Haimanll2003l;lRaskutti et alj|2012l) . 

In the simplest picture, the competition between photo-heating 
and cooling due to the adiabatic expansion of the Universe results 
in a power-law temperature-density relation following reionisation, 
T = ToA''"^ for A = p/{p) < 10. The density dependence of the 
recombination rate means that higher density gas recombines faster. 
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yielding more neutral atoms per unit time for photo-heating. These 
regions thus cool less rapidly than lower density gas, resulting 
in a temperature-density relation which evolves from isothermal 
(7 ~ 1) following re ionisation toward a power-law with 7 ~ 1.6 
l Hui&Gnedirilll997h . n principle, if the redshift dependence of 
the temperature-density relation can be measured then the timing 
of reionisation can thus be inferred. Indeed, by stud ying the ther- 
mal w idths of absorption lines in the Lyman-Q forest. ISchave et aU 
( I2OOQ) observed an increase in the temperature at mean density and 
a flattening of the temperature-density relation at z ~ 3, which may 
indicate the epoch of He 11 reionisation occurred around this time. 

However, this is a somewhat simpli fied picture; accordin g 
to the numerical simulations presented bv iMcOuirm et al. 
He II reionisation is inhomogeneous and long-range heating by 
hard photons will induce large-scale fluctuations of the order of 50 
comoving Mpc in the IGM temperature and He 11 ionisation state. 
Another question mark hangs over the slope of the temperature- 
density rel ation describ i ng the IGM therma l state. Obs ervational 
work f rom lBecker et al.l ( l2007h , lBolton et alj jZOOSi) and lViel et all 
( l2009h has suggested that the IGM may obey an 'inverted' 
(7 < 1) temperature-density relation in which, somewhat counter- 
intuitively, less dense gas is hotter than denser gas. Although it 
appears difficult to produce this result by H e 11 photo-heating by 



quasars jMcOuirm et al. 20091: 



been suggested bv lChang et al 



Bolton et an[200a). it has recently 
( I2OI ih and lPuchwein et al.l ^201 11) 



that it could be a consequence of volumetric heating by TeV emis- 
sion from blazars. 

Given these uncertainties, it is important to investigate the ob- 
servational constraints in more detail. One way to achieve this is 
to directly compare different methods used to measure the IGM 
thermal state. This allows one to establish whether these different 
approaches are consistent, or whether there are systematic uncer- 
tainties which impact differently upon the competing approaches. 
The methodologies used in the literature thus far to measure the 
IGM thermal state from the Lyman-a forest can be broadly divided 
in two classes. The common feature in both approaches is that they 
measure the IGM temperature via the impact of Jeans (pressure) 
smoothing and thermal Doppler broadening on the Lyman-a for- 
est. The first class consists of methods that fit Voigt profiles to 
each absorption line in t he Lyman-g f orest. Example s of th i s class 
are the earlier work of ISchave etalj |2000), Ricot ti et alJ j2000h 
and lMcDonald et alJfeOOOl) and more recentlv Bolto n et alTj20IOl) 
The second class consists of methods in which the transmitted flux 
is analysed with a global statistical approach, without decompos- 
ing the spectr a into separate features. Power spectra stu dies belong 
to this class fealdarriaga et alj|200ll ; fviel et alj|2009h . as well as 
methods that examine other statistical properties of th e forest, such 
as th e flux probability d istribution function (PDF) jBolton et al.l 
2008; Caluraet al. '2012*), the wavelet analysis meth od applied in 
Lidz et al., ( i2010i) and the 'curvature' statistic used bv lBecker et al.l 



( 120 III) . 

The main aim of this work is to compare two of these com- 
peting techniques, the flux PDF and wavelets, by applyi ng them to 
the m etal-cleaned Lyman-a forest spectra presented bv lKim et al.l 
( l2007i) . These data have previously been used in studies of the 
flux PDF which have found the IGM temperature- density rela- 
tion may be isothermal or inverted at 2 < 2 < 3 jBolton et al.l 
[2OO8; Viel et al. 2009). In order to interpret the results we have also 
utilised and exten ded the suite of hyd rodynamical simulations used 
in the analysis of iBecker et alj 1 I2OI li) . T his compari s on is of par- 
ticular interest because, as pointed out bv lLidz et al.l j20IOl) . there 
appears to be some tension between the IGM thermal parameters 



inferred from wavelets and the flux PDF, particularly with respect 
to measurements of the slope of the temperature-density relation. 
This paper is therefore organised as follows: in Section 2 we re- 
view the observations and numerical simulations used; in Section 
3 our implementation of the wavelet analysis is presented; in Sec- 
tion 4 we describe our interpolation and parameter determination 
methodology and in Section 5 we present our results along with a 
comparison to previous studies. We conclude in Section 6. An ap- 
pendix at the end of the paper lists the simulations we use in this 
study, along with several tests for numerical convergence and the 
sensitivity of our results to parameter assumptions. 



2 OBSERVATIONS AND NUMERICAL SIMULATIONS 

In this section we briefly review our observational dataset, the nu- 
merical simulations used for our theoretical interpretation and the 
procedure used to generate synthetic spectra. 



2.1 Spectra, metal removal and continuum placement 

Our an alysis uses the eighteen metal-cleaned quasar spectra from 
iKim e t al. (2007) obtained with the Ultraviolet and Visual Echelle 
Spectrograph (UVES) on the VLT. The spectra are sampled with 
pixels of width 0.05 A and have a signal to noise ratio per pixel 
of order 30-50. In order to avoid the proximity effect, the region 
4000 kms~^ bluewards of the Lyman-a emission line has been 
excluded. These spectra contain Lyman limit systems with column 
densities in the interval 10^^'^cm^^ ^ A'^hi ^ lO^^cm^"^, but no 
damped Lyman-a systems, defined by a column density Nm ^ 
lO^O'^cm-^. 

In order to correct the H i absorption for metal contamination, 
metal absorption lines in the spectra were identified and fitted with 
Voigt profiles. These were substituted by a continuum level or by a 



Lyman-a only abs orption profile g enerated from the fitted Lyman- 
a parameters (see lKim et alj|2007l for details). Note this approac h 
differs from the metal removal procedure used in lLidz et al.l ( l201(ih . 
where only narrow absorption lines (with b < 7km/s) were iden- 
tified as metals and excised from the spectra. 

The continuum level in the spectra was determined by locally 
connecting regions that are thought to be absorption-free. This is 
an iterative procedure, which starts with connecting non absorbed 
regions and is subsequently updated during the process of Voigt 
profile fitting the H i and metal lines. Note that in our analysis we 
neglect the possibility of having an extended and slowly varying 
continuum absorption. This means that spectral regions which are 
considered to be absorption free could actually suffer from absorp- 
tion by a broad H i density fluctuation, and the measured optical 
depth would then be underestimated. We shall discuss the effect of 
continuum placement on our results further in Section lSTI 



2.2 Numerical simulations 

T he simulations we u se are based on the suite of models used 
bv lBecker etai] ( l20Ilh . which we have extended by further varying 
the cosmological and IGM thermal parameters assumed. We make 
use of Gadget-3 code w hich is a parallel smoothed particle hy- 
drodynamics code dSpringe l 2005). Our simulations are performed 
in a periodic box of 10 comoving Mpc/h in linear size. We de- 
scribe the evolution of both the dark matter and the gas, using 256"^ 
particles for simulations in which the cosmological parameters are 
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Table 1. Properties of the quasar spectra from lKim et alj J2007h . Zcm is the 
approximate redshift of the quasar measured from the Lyman-a emission 
line; ^Lya and \\^yoi are the redshift and wavelength intervals associated 
with the Lyman-o absoiption; S/N is the signal to noise ratio per 0.05A 
pixel. 



varied, or 2 x 512'^ particles for simulations in which the IGM ther- 
mal state parameters are varied. A summary of the simulations is 
given in Tab lAll in the appendix. 

The simul ations all start at z = 99 with initial conditions gen- 
erated using the lEisenstein & Hul ( Il999l) transfer function. Star for- 
mation is incorporated using a simplified prescription in which all 
gas particles with A > 10^ and temperature T < 10^ K are con- 
verted into coUisionless stars. Since we are not interested in the de- 
tails of star formation, we thereby avoid the small dynamical times 
that would arise due to these overdense regions. As the bulk of the 
Lyman-Q forest absorption corresponds to densities A < 10, this 
prescription has an impact at below the percent level on the final 
computation of fiux PDF and flux power t Viel et al. 2004 ) and so 
we expect this will not affect our work. To check numerical con- 
vergence of our simulations we have also performed a series of 
simulations with varying gas particle and box size. We conclude 
from these tests, demonstrated in Figure lA2l in the appendix, that 
our study of the statistics of small-scale structure of the Lyman- 
Q forest demand the relatively high mass resolution afforded by a 
large number of particles, 2 x 512^, in a 10^^^ Mpc volume. With 
greater computational resources, we could improve the accuracy of 
our simulations by increasing the box size from 10 h^^ Mpc, while 
maintaining high mass resolution. However we believe that the pos- 
sible improvements are small compared to the error budget we have 
assumed in Section [J!2l 

A spatially uniform ultraviolet (UV) background applied in 
the optically thin limit determines the photo-heating and photo- 
ionisation of the gas in the simulations. In order to generate dif- 
ferent thermal histories, we have re scaled all of the H i , He i and 
He II photo-heating rates using the 'Haardt & MadatJ j200lh UV 
background model for galaxy and QSO emission. The photo- 
heating rates, ti, were changed by rescaling their values in a density 
dependent fashion, ef°" = CA«ep'^°^ A list of the values used 
for the scaling coefficients and 5, along with the correspond- 
ing values for To and 7 and the other main simulations parame- 
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Figure 1. Redshift range of the observed spectra listed in Table 1 . The two 
vertical lines indicate the edges of the redshift bins with effective redshifts 
of (2) = [2.07, 2.52, 2.93], and which are roughly centred on the redshift 
of our simulation outputs. The percentage of spectra in each redshift bin are 
57, 25 and 18 per cent, respectively. 



ters are given in Table IaTI in the appendix. The resulting thermal 
histories are self-consistent in the sense that the gas pressure, and 
hence the Jeans smoothing, is compatible with the gas temperature. 
Computing simulated spectra with this treatment requires a large 
computational effort, because e ach thermal history requires its own 
simulation. By comparison, the lLidz et al.l (2010) reference simu- 
lation was run with a fixed chosen ionisation state, and the thermal 
state was then superimposed by applying a temperature-density re- 
lation in a post-processing step, thereby neglecting a full treatment 
of Jeans smoothing in this approximation. However, note that both 
our approach and the Lidz et al. (2010) technique ne glect radiative 
transfer effects (e.g. Titfle v & Meiksinll2007i : iMcOuinn et al.l2009h 
because of the computational effort that solving the cosmological 
radiative transfer equation would imply. 

We have selected three simulation snapshots at 2; = 
[2.17,2.55,2.98] from the Becker et al. (2011) models, which 
cover the redshift range of our data, and the values of To and 7 
have been determined by fitting the temperature-density relation 
at each redshift for each simulation. These snapshots then deter- 
mine the way in which we split the data into three redshift ranges, 
divided at z = [2.35,2.70], as shown in Figure [T] These three 
redshift bins contain 57, 25 and 18 percent of the data, respec- 
tively, and the effective (average) redshift of the data in each bin is 
(z) = [2.07, 2.52, 2.93]. Hence-forward, these effective redshifts 
will be the nominal values used when quoting our results. 

In order to assess the effect of astrophysical and cosmologi- 
cal uncertainties on the IGM physics we vary the simulation pa- 
rameters on a grid, one at a time. There are two sets of simu- 
lations: in the first set we vary only the IGM thermal state pa- 
rameters around a reference model with [f2ni, f^b, ^, o"8, ?isl = 
[0.26, 0.0444, 0.72, 0.8, 0.96] and with a reference IGM thermal 
state [ro/10^K,7] = [15,1.6]. The range of thermal state pa- 
rameters covered by our simulations extend from around To = 
4 600 - 31 000 K and 7 = 0.7 - 1.6. In the second set we 
vary only the cosmological parameters around [fim, ft, cs, ns] = 
[0.26, 0.72, 0.85, 0.95] ± [0.04, 0.08, 0.05, 0.05], with a reference 
IGM thermal state [To/10^K,7] = [20,1.06] and fixed Sib = 
0.0444. The difference in cosmological parameters for these two 
simulations sets owes to imperfect planning, though this has not 
been a great source of concern or bias since, as we demonstrate in 
the appendix of this paper, the uncertainties in the cosmological pa- 
rameters are not the limiting factor in our predictions for Lyman-a 
spectra or in our interpretation of the data. 
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Figure 2. Real part of the Morlet wavelet. Equation (ij, with smoothing 
scale s = 70kms~^. The width of the envelope depends on the scale 
probed by the oscillations. 



2.3 Synthetic spectra 

Our synthetic spectra are obtained usin g the following proc edure 
described in iTheuns et alJ h998h and lBoltonet alJ j2008h . and 
briefly reviewed here. At each redshift slice, approximately 10'' 
randomly chosen lines of sight are selected. We convolve the re- 
sulting H I de nsity with a V oigt profile using the approximation 
introduced in IXepper-Garcfal yOOfi) - this approximation is suffi- 
ciently accurate for the range of column densities considered here. 
We resample the spectra into velocity intervals of 4.4kms~^ and 
we account for instrumental resolution by convolving the spectra 
with a Gaussian with full width at half maximum of 7 km s^^ . We 
add noise to the spectra with the same level of the observed data. 

In order to leave the Hi effective optical depth, Tcb, as a free 
parameter we adjust the mean transmitted flux as (F) — > {F') = 
(exp [— T^]), where A is chosen such that the global mean normal- 
ized flux {F') satisfies Tcfi = — ln((F')). As we will discuss, Tcft 
is one of the main uncertainties limiting the determination of the 
IGM thermal state in our wavelet analysis. 



3 WAVELET ANALYSIS OF SPECTRA 

The wav elet decomposition of the Lyman-g forest was first pro- 
posed bylTheuns & Zaroubi' fe OOOh . lMeiksinI JlOOOh . lTheuns et"ai] 
( l2002n and Zaldarriaga (2002). The idea is to filter the spectra 
in such a way as to construct observables that are sensitive to 
the thermal state of the IGM. The wavelet decomposition was 
therefore suggested because it might be used to detect temporal 
or ev en spatial variations of physical properties of the IGM (e.g. 
iMcOu inn et al. 2011). 

The physical motivation for this type of filtering relates to the 
effects of Doppler broadening and Jeans smoothing on small scale 
structure in the Lyman-a forest. The thermal Doppler effect arises 
due to the velocity dispersion of the hydrogen atoms, which causes 
broadening of the absorption lines ( f or a recent analysis of thermal 
broadening and Jeans smoothing, see lPeeples et al. l( l2010allbl) ).The 
velocity distribution is described by the Maxwell distribution 



P{v)dv - 



2TTk-BT 



exp 



2k^T 



dv. 



(1) 



from which it can be seen that the velocity dispersion is propor- 
tional to Vt. This results in a broadening of the absorption spectra 
by a factor y^^^ ^ 13 km s"^ at T = 10'' K. The other impor- 
tant effect is Jeans smoothing. Because the ideal equation of state 
holds, an increase in the temperature of the gas corresponds to an 



increase in pressure w hich then also s mooth es the small-scale struc- 
tures (see for example [Pawlik et al.l | |2009|) 1. The Jeans smoothing 
depends on the full thermal history of the IGM, because pressure 
forces alter the dynamical state of the gas fflui & Gnedia,1997.) . 



3.1 The wavelet amplitude PDF 

Following iLidz et ah I ( l2010h . we have implemented the 'Morlet 
wavelet' as a probe of this thermal smoothing, which we briefly re- 
view here. The Morlet wavelet is a Gaussian in the complex plane, 
defined in velocity space by 



i^k{v) = .4 exp [—ikv/2TT] exp [— u^/2s'^] 



(2) 



where the 'smoothing scale' s — 27r/fc is chosen so that the 
wavelet changes its global width depending on the probed scale 
k. Our normalization constant A is chosen so that the integral of 
the squared wavelet function is unity. An example is shown in Fig- 
ur e[2]for a smoothin g scale s = 70 km s~ ^ . This follows the choice 
of iLidz et ai]( l2010() which taken as a compromise between max- 
imising the sensitivity to the small-scale structure and avoiding the 
possible contamination by metal lines. 

The spectra F{v) are first convolved with the Morlet wavelet, 
to obtain a filtered spectrum 



f{v)= / dv' F{v)^p{v — v'). 



(3) 



The filtered spectrum / is then squared and smoothed in order to 
compute the 'wavelet amplitude' 



(4) 



A(«) = i I dv'e{\v-v'\,L/2)f{v'), 

where 0{v,L/2) is the top-hat function with width 
L — lOOOkins^^. This choice of large-scale smoothing 
follows iLidz et alj bOld) , though we have also checked that our 
results do not depend strongly on the exact value of L. 

Figure [3] shows an example of these processing steps for 
Q0055-269, our highest redshift spectrum. It can be seen that the 
wavelet amplitude is greater in regions of the spectrum with ab- 
sorption lines of width comparable to the sampling scale. The re- 
sulting signal captures some of the inhomogeneity of the original 
spectrum. 

The final observable is the 'wavelet amplitude PDF', p{A), 
which is calculated by binning A, calculated for each spectra, into 
a histogram with ten logarithmically spaced bins, whose minimum 
and maximum values are taken from the data. The same processing 
is applied to our simulations to obtain the predicted wavelet ampli- 
tude PDF. The resulting wavelet amplitude PDFs computed from 
our data are presented in Figure |4] These data are the key observa- 
tional results of this work, which will be analysed and interpreted 
in Section|5] 



3.2 Error bar estimates 

Following iLidz et ah I 1 I2OO6II we have attempted to estimate the 
error bars of the wavelet amplitude PDF using the 'jackknife' 
method, although as we describe below, we believe our main un- 
certainty arises from the accuracy with which we can predict the 
wavelet amplitude PDF. The jackknife method is a resampling 
method in which our spectra are divided into ng = 10 subgroups of 
equal size, from which a set of rig wavelet PDFs, pk{A), are com- 
puted by omitting from the data one subgroup of data at a time. The 
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Figure 3. Upper panel: the original QSO spectrum, Q0055-269. Middle panel: the amplitude of the signal, Equation ^3)> obtained by convolving the spectrum 
with a Morlet wavelet with s = 70kms~^. Lower panel: the same amplitude after squaring and smoothing with a lOOOkms"^ wide top-hat filter. 
Equation (4). 
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Figure 4. The observed wavelet amplitude PDF. Equation @, for s = 70kms^^ and L = 1 000 kms^^, for redshifts z = [2.1, 2.5, 2.9] (left to right). 
The error bars displayed here, calculated via jackknife resampHng with Equation (5), are thought to be underestimated and so we have made an allowance for 
theoretical uncertainty in our interpretation of the data, as described in Section [331 



covariance matrix of the wavelet PDF amplitude is then estimated 
using 

rig 

= ^[piA,) - MA^)MA,) ^ p,{A,)]. (5) 
fc=i 

iLidz et all ( l2006h tested their covariance values against those es- 
timated directly from 10 000 mock spectra, and found that Equa- 
tion ^ holds approximately, but that the error bars can sometime 
be underestimated, especially in the tails of the wavelet amplitude 
distribution. Partly owing to this observation, and partly due to our 
more limited number of mock QSO spectra on which to perform 
the jackknife method, we have opted for what should be a con- 
servative estimate of wavelet PDF uncertainties: the diagonal ele- 
ments of Equation ([5j are all replaced with (0.25 x ma,x{p{Ai))^ , 



as shown in Figures [6] and [T] and we ignore the off-diagonal ele- 
ments suggested by the jackknife resampling. In doing so we are 
attempting to make an allowance for the theoretical uncertainty as- 
sociated with our calculation of the wavelet amplitude PDF as well 
as for the level of interpolation errors suggested by our validation 
tests described in Section \4J\ We note that we implemented the 
jackknife method, as defined in Equation (O, (as well as a boot- 
strap resampling method) to estimate the covariance and error bars 
shown in Figure|4l but we concluded that the covariance was being 
underestimated. Given the approximation applied in our interpola- 
tion scheme, we believe the current error-bar prescription we ap- 
ply is sufficiently representative to constrain the central value and 
width of the wavelet amplitude PDF, and that it is on the conserva- 
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live side. A full theoretical understanding of the wavelet amplitude 
bin-bin covariance remains an open problem. 



4 PARAMETER DETERMINATION METHODOLOGY 

In this section we describe the core ingredients of our analysis: 
our interpolation scheme, parameter sampling method and IGM 
parametrization. 



4.1 Interpolation scheme 

In order to calculate the wavelet PDF for a given location in our 
parameter space, we perform an interpolation of the wavelet PDF 
calculated over our available simulations. Our approach is to per- 
form a cubic-spline interpolation of the wavelet PDF differences 
as a function of the cosmological and astrophysical parameters, for 
which we have three simulations along each cosmological parame- 
ter direction, three simulations varying 7 and six simulations vary- 
ing To, as illustrated in Figure[5] 

As mentioned in Section [23] the effect of varying Tcfi is cal- 
culated in post-processing. We therefore calculate the wavelet PDF 
on a fixed grid of 100 values - checking that an implementation 
with ten grid points in this direction would also be acceptable. We 
then apply the scheme 



p(A){eo,T^a) + 

n 



,do},T^a), (6) 



for the two closest values of Tea on our grid and obtain the final 
p{A){9,Tch) via linear interpolation of these two function eval- 
uations. Here 6 denotes a parameter vector with n components, 
and Ap is the interpolated PDF differences relative to the fiducial 
model, 60. 

We have checked the accuracy of this interpolation scheme by 
comparing it with the wavelet PDF calculated for a 'validation set' 
of simulations which are not used in the interpolation procedure. 
The wavelet PDF for these simulations, ('DIO' and 'ElO' with pa- 
rameters given in Table I A it were found to be satisfactorily repro- 
duced, within the uncertainties that we have assumed, as described 
in Section \32\ An illustrative case is shown in Figure [6] for the 
z = 2.9 bin. 

We believe that this interpolation captures the variations in the 
wavelet amplitude PDF predicted by our simulations sufficiently 
accurately for our purposes and within the generous errors we have 
assumed. We therefore expect that our parameter constraints will be 
on conservative side. We suggest that, with more computational re- 
sources than currently available to us, the PICO jFendt & Wandelj 
[2007) training-set/interpolation scheme could be applied to accu- 
rately and efficiently predict Lyman-Q forest observational quanti- 
ties. 



4.2 Parameter sampling 

In order to estimate the astrophysical and cosmological parameters 
and their uncertainties, we use a sampling-based approach. Bayes' 
theorem can be written as £{'D\6)n{9) — ZP(6), where £ is the 
likelihood of the data D given parameters 6, ti(9) are priors on the 
parameters, Z = J d6n{9)C is the model likelihood or 'evidence' 
and P(9) is the sought-after posterior distribution of the parame- 
ters. 
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Figure 5. Simulations used in our interpolation scheme: reference model 
(crosses), Tq interpolation (circles), 7 inteipolation (squares) and the val- 
idation set (diamonds). Further information about the simulations may be 
found in Table lAl] 



We perform posterior sampling using the 'nested sain pling' 
algorithm, which is technique proposed by ISkillind Jiooi) . prin- 
cipally for estimation of the evidence Z with posterior sampling 
as a by-product. We have chosen this algorithm because it is well- 
suited for sampling likelihoods that are multi-modal, strongly de- 
generate, or non-Gaussian. To briefly summarise the nested sam- 
pling method, the multi-dimensional integral Z is remapped onto a 
particular one dimensional integral. A key ingredient for this kind 
of sampling is then the ability to draw uniform random samples 
from the remaining region of parameter space delimited by an iso- 
likelihood surface. We made use of t he publicl y available imple - 
mentation MultiNest jperoz & Hobs on 2008; Feroz et al.l2009t) . 
in which those regions bounded by iso-like lihood contours are ap- 
proximated by ellipsoids (Mukheriee et aljr2006,) . We run MULTI- 
Nest in a configuration with 500-1000 'live points' and with a rel- 
atively low sampling efficiency of around 10"''. The posterior sam- 
ples are then analysed with the GETDIST package (L ewis & Bridie 
I2OO2I) in order to extract two-dimensional and one-dimensional 
marginalised constraints. 



4.3 IGM parametrization and priors 

The next step in our analysis is to parametrize the redshift evolu- 
tion of the quantities we wish to constrain. We therefore investigate 
a model for the possible redshift dependence of the IGM parame- 
ters in which To, 7 and are allowed to vary as piecewise con- 
stants centred on the redshifts z = [2.1, 2.5, 2.9], and which will 
be referred to as the 'redshift bins parametrization'. 

We have put wide flat priors on the IGM parameter ranges, 
5 < To/IO^K < 50, 0.5 < 7 < 1.7, 0.075 < r^alz = 2.1] < 
0.2, 0.13 < Tes[z = 2.5] < 0.3, and 0.2 < t^sIz = 2.9] < 0.5; 
here the Tofi prior appro ximates and enco mpasses the ranges al- 
lowed from Figure 13 of iKim et al.l ( |2007|) . For the cosmological 
parameters, we have imposed the flat priors 0.6 < ag < 1.0, 
0.9 < Us < 1 and 0.20 < < 0.32, and 60 < -ffo < 
84 which is intended to be a conservative range encompassing 
the region favoured by WMA P jLarson et al]|201 ll) and Ho con- 
straints jFreedman etai]l200ll) . 

Note that more restrictive power-law parametrizations for 
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sorption at 2: < 3 (e.g. iBecker et aljbol ih . Increasing Tch lowers 
the characteristic gas density probed by the Lyman-Q absorption, 
shifting the peak amplitude of the wavelet PDF to higher values 
due to the colder underdense gas present for our fiducial 7 = 1.6. 
Clearly the uncertainty in Tcfi needs to be marginalised over in or- 
der to estimate To and 7, as any physical effect that affects the 
power spectrum of the Lyman-Q absorption lines at small scales 
fc > 0.1 s/km will also substantially affect the wavelet amplitude 
PDF. 

Finally, the effect of the cosmological parameters on the 
wavelet amplitude is found to be weak, except for erg that has an 
slight impact on the wavelet PDF at redshift 2; = 2.1 as we demon- 
strate in Figure IaTI of Appendix lAl I The effect appears not to be 
as simple as a shift in the wavelet amplitude PDF as m ight be ex- 
pected for a change in the power spectrum, as argued bv lLidz et al.l 
bold) . Our explicit simulations of the effect of varying erg show a 
change in the width of the wavelet PDF in the lower redshift bin, 
where structure formation is more advanced. 
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Figure 6. Comparison of the wavelet PDF calculated with our interpola- 
tion scheme (dotted) to the PDF from a direct simulation (solid), for two 
validation simulations, CIO (upper) and ElO (lower). The agreement is sat- 
isfactorily within the uncertainties that we have assumed. 



5 RESULTS 

We now turn to describing the results of our analysis. Briefly, to 
restate the main aim of this work, we wish to compare the IGM 
cosmological and astrophysical const raints derived fro m our re- 
analysis of the Lyman-a flux PDF of iKim et al.l ( 1200 7') with our 
new constraints derived from the wavelet amplitude PDF. We ap- 
ply these two methodologies to the same dataset for the first time, 
enabling us to explore any systematic differences between these ap- 
proaches. The com putation of the flu x PDF likelihood is performed 
in a similar way to lViel et al.l ( |2009|) . using a second-order Taylor 
expansion of the cosm ological and astro physical param eter space. 
We have extended the IViel et al.l JioO^ analysis of the iKim et al.l 
( l2007h PDF to use our redshift bins parametrization, thereby di- 
viding the data set into three redshift bins as opposed t o using a 
sin gle powe r-law parametrization as in IViel et alj (I2OO?). Follow- 



ing [\2^1_etal. (2009), we also fit the simulations to the observed 
flux PDF in the flux range F = 0.1-0.8 only, in order to minimise 
the effect of continuum uncertainties on our results. 



To(z) , 7(2) and Tcsjz) have been inves tigated by a number of au- 
thors dViel et al.ll2009l: iBird et alJlioTH) . Owing to the potentially 
complex and uncertain phenomenology suggested by theoretical 
models of the IGM, we have opted for the more general redshift 
bins parametrization. 

4.4 The effect of astrophysical and cosmological parameters 
on the wavelet PDF 

Having described our IGM parametrization and implemented our 
interpolation scheme, we are now in a position to demonstrate how 
the wavelet amplitude PDF depends on the astrophysical param- 
eters. To, 7 and t^s, before proceeding to present our tempera- 
ture constraints. This is illustrated in Figure|7] from which we may 
co nfirm the phenom enology for the wavelet amplitude PDF found 
m iLidz et alj jloToh : the wavelet amplitude PDF shifts to smaller 
values for higher temperatures owing to Doppler broadening and 
Jeans smoothing which suppresses small scale power. Higher val- 
ues of 7 also shift the peak of the wavelet PDF to lower amplitudes 
at fixed To due to the decreased thermal smoothing associated with 
absorption from overdense gas, which dominates the Lyman-Q ab- 



5.1 Constraints on the IGM thermal history 

Our main results are shown in Figure [8] which displays the Icr and 
2(7 contours obtained from our analysis of the wavelet PDF (yel- 
low contours) and the flux PDF (red contours). Firstly, we note that 
the values of To and 7 inferred from our analysis of the wavelet 
and flux PDF are in broad agreement with each other, and are for- 
mally consistent within la. In part this is because in this analysis 
we have, in the first instance, left t^s as free parameter, which en- 
larges the parameter space consistent with the wavelet PDF due to 
the degeneracy between t^s and To. It is also apparent, however, 
that the flux PDF constraints generally favour a lower value for the 
temperature-densi ty relation slope, 7, co mpared to the wavelet PDF 
(e.g. lBolton et aLl2008l : IViel et al.l2009l) . This implies there may be 
a systematic difference between the two methodologies. 

Motivated by recent suggestions about possibl e systemat- 
ics in the flux PDF due to continuum placement dLed l20Ilt 
iMcOuinn et al.l201 1)) we have investigated whether the wavelet and 
flux PDF show any sensitivity to the continuum level assumed. We 
performed a check by low ering the continuum level on the synthetic 
spect ra by 3 per cent (e.g. lTvtler et alj2004l : lFaucher-Giguere et al.l 
|2008|) . and then recalculating both the flux and wavelet PDF for 
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Figure 7. The predicted wavelet amplitude PDF at redshifts z = [2.1, 2.5, 2.9] (left to right columns) for our fiducial cosmological parameters. Top panels: 
varying Tq. The wavelet PDF shifts to lower values for higher temperatures owing to the suppression of power at small scales. Middle panels: varying 7. 
Here 7 defines the slope of the temperature-density relation, and a higher value for 7 leads to a higher value for the temperature at fixed Tq in the overdense 
gas predominantly probed by the Lyman-a forest at 2 < 3. Bottom panels: varying T^g. Increasing T^ff decreases the mean transmitted flux and shifts the 
characteristic density probed by the Lyman-o forest to lower values. These regions are cooler for our fiducial 7 = 1.6 than more overdense regions, leading 
to shift towards higher values of the wavelet amplitude. 



our 2 = 2.98 simulations, for two models with ToflO^K] = 18, 
Teft = 0.350 and 7 — 1.6 and 0.7 respectively. Figures showing 
the results from these tests are shown in Figure lASl in the appendix. 
We conclude that while the flux PDF shows some sensitivity to the 
continuum level at high (F > 0.8) and low {F < 0.1) fluxes (data 
that have not been used in our fit), the overall change in the fit is 
small compared to the improvement in the fit from moving (for 
example) from 7 — 1.6 to 7 = 0.7. It therefore appears that con- 
tinuum errors do not fully explain the tendency for the flux PDF to 
favour somewhat lower values of 7 at z = 3. 

For comparison, we also show i n Figu re [S] results from the 
wavelet PDF analysis of Lidz et all 1 2010l) extracted from their 



Figure 25. The three iLidz et alj ilOKf) redshift bins are z = 
[2.1, 2.6, 3.0] which rough ly correspond to our own redshift bins. 
Note that iLidzet all ilOldl imposed a prior 7 — 1.0—1.6 and 
regarded their constraints as approximate, cautioning against tak- 
ing their results too literally. Therefore, in interpreting their results 
we will attempt to make an allowance for the caveat they expressed 



by examining their 2a constraints. iLidz et al.l feOlOh have analysed 
just over double the numb er of spectra used in thi s study by using 
the 40 spectra reduced bv iDaH'Aglio et all feOOSi) : sixteen of our 
eighteen spectra also appear in their sample. Our wavelet PDF er- 
ror bar treatment is probably the more conservative, explaining why 
our constraints are somewhat looser. In genera l, however, we find 
very good agreement with the lLidz et alj ( l201(il) constraints, which 
is encouraging given the independently reduced observational data 
set and different simulation method used by these authors. 

Our one dimensional parameter constraints are shown in Fig- 
ure |9| in which we can compare the results from the wavelet and 
flux PDF side by side. It is clear that the cosmological parame- 
ters are unconstrained by the wavelet PDF owing to their weak 
effect, with perhaps only a lower bound on as being found. The 
flux PDF - a one point statistic of the unfiltered spectra - puts 
the stronger constraint on the values of Tbs; the values of T^ff 
we find for the se redshift bins agree at the l -2a level with those 
determined by iFaucher-Giguere et alj ( |2008[) . This is in fact our 
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Figure 8. Upper panels: constraints on the IGM the rmal parameters T and 7 from our implementation of wavelet PDF analysis (yell ow, light fil l ed con tours), 
compared with our reanalysis of the flux PDF from lViel et al.l )2009l) (red, dark filled contours) and the Icr and 2(7 constraints from lLidz et al.l )201(j) which 
assume a prior -y = 1.0 — 1.6 (black contours). Lower panels: the joint constraints from the wavelet and flux PDF (filled contours) again compared to the 
iLidz et all feoiOl) ' 



I wavelet PDF constraints. 



main motivation for performing a 'joint analysis' of the flux PDF 
and wavelet PDF together (combining their Ukehhoods with equal 
weight), in order to self-consistently add the Tcff constraint derived 
from the flux PDF onto the parameter space consistent with the 
observed wavelet PDF. The final results for the joint analysis are 
summarised in Table |2] and the lower three panels of Figure [8] The 
joint constraints favour an IGM temperature-density relation which 
is close to isothermal, with temperatures at mean density which lie 
in the range 10 000-20 000 K. However, the la uncertainties are 
too large to infer any significant redshift evolution in these quanti- 
ties. 



5.2 Comparison with previous constraints 

Finally, we compare our results to previous constraints in the lit- 
erature and consider the implications for the thermal history of the 
IGM. In Figure[TO]we present a comparison of our results to models 
for the redshift evolution of the temperature at mean density, To(z), 
from the literature. We show our joint constraints on To (black error 
bars) and from the wavelet PDF alone (red error bars) t o gether with 
the models of Hen reionisation from McOuinn et al. (20091) (up- 
per panel) and blazar heating from Puchwein et al. c201 ih (lower 
panel). Specifically, we have c ompared the available data with 
the 'Dr and 'S4b' models from lMcOuirm et ( l2009h : the latter 
model implements a harder quasar UV spectral index of 0.6 com- 
pared to their fiducial 1.6. Our constraints are generally consistent 
with the iMcOuinn et alj ( l2009t) model of He 11 reionisation, with 
the softer UV spectral index model preferred. 

The wavelet only constraints tend to favour the weak blazar 
heating model of IPuchwein et al.l ( 1201 ll) . whose median tempera- 
ture at mean density is shown. In contrast, IPuchwein et al.l j201ll) 
concluded that their intermediate blazar heating was preferred 
based on a comparison of the temperature of their models calcu- 
lated at the 'optimal overdensity' (which rises to A ~ 6 at 2: = 2) 



probed by the ISecker et al.l ( l20Ilh curvature constraints. The dif- 
ferences between the blazar heating models are more pronounced 
at mean density ; the temperature-density relations predicted by 
Puch wein et al.l (2011) are simil ar for all model s at A ^ 2-3, 
which partially accounts for why IPuchwein et al] j201 ih conclude 
the intermediate model is favoured. We note, however, the wavelet 
PDF is sensitive to gas temperatures over a range of densities (in- 
cluding the mean density) as it is a distribution rather than a single 
number (i.e. the curvature statistic used by Bec ker et al. 201 1). A 
precise measurement of the temperature-density relation could in 
principle rule out the blazar heating model if 7 > 1, as the vol- 
umetric heating rate used in the blazar heating models produce a 
strongly inverted temperature-density relation hy z = 2. 

We also compa re our constraints on To with the measurements 
of lLidzetJ] ( l2010l) . also u sing the wavelet PDF technique, and the 
curvature measurements of ISecker et alj ( l201lh (plotted assuming 
7 = 1.3, which is within ~ 1.5a of our joint constraints in all 
redshift bins). Note that although they used a different measure- 
ment technique and data, ^ Becker et al.l ( I2OI 1*) use the same set of 
hydrodynamical simulations as us. There is good agreement be- 
tween the measurements at z < 3, although there doe s appear to 
be som e tension between the iBecker et alj ( l201lh and iLidz et all 
( I2OIOI) results in the r edshift range = 3 — 3.5. As pointed out 
bv lBeckere"tal] ( l201ll) . however, differences in the effective opti- 
cal depth assumed in the two studies may play a role here; as we 
have demonstrate d the wavelet PDF is rather sensitive to Tcff. The 
iBecker et al. Ij) constraints also have significantly smaller error 
bar s compared to the w avelet PDF measurements. Note, however, 
the IBecker et al.l ( 12011 *) measurements do not attempt to simulta- 
neously measure both To and 7. As noted previously, they instead 
measure the IGM temperature at the characteristic density, A{z), 
probed by the Lyman-Q forest. Their constraints do not marginalise 
over the uncertain values of 7 and Tofi, and so a direct comparison 
of their uncertainties to our results is less straightforward. 
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at z S> 3. Overall, we find our results are in agreement with 
previous studies of the IGM temperature that suggest there may 
be additional heating in the IGM at z < 4, most likely due to 
the reionisation of He ii b y quasars (see also recent studies of the 
He II Lyman-g fore st, e.g. IShuU et allllOlol : IWorseck et allllOllI : 
ISvphers et al.ll20ll]) . We conclude improved measurements of the 
slope of the temperature-density relation will be required for testing 
blazar heating models in detail. 



6 CONCLUSIONS 

We have analyzed eighteen metal cleaned Lyman-a forest spectra 
in the redshift range 1.7-3.2 using both the flux PDF and wavelet 
PDF. The results have been interpreted using a suite of hydrody- 
namical simulations to place constraints on the thermal state of the 
intergalactic medium while marginalizing over the uncertainty in 
the cosmological p arameters. Our wav elet analysis is similar to the 
analysis performed'Xheuns et al. ( 2002), following most closely the 
technique employed by Lidz et al. (2010), but using independently 
reduced spectra and a different approach to simulating the IGM. An 
analysis of the constraints on the IGM thermal state obtained from 
the flux PDF using the same data set furthermore enables us to ex- 
plore any systematic differences between the two methodologies. 
The main results of our study are as follows: 

• The constraints on the IGM thermal state al z — 1.7-3.2 
derived from the wavelet PDF and flux PDF analysis are formally 
consistent with each other within the rather large uncertainties. 
However, we find there is some mild tension between the two mea- 
surements, with the flux PDF measurements generally preferring a 
lower value for the slope of the temperature-density relation at all 
redshifts. 
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Figure 9. One dimensional marginalised constraints, coiresponding to the 
two dimensional constraints shown in Figure [8] Note the constraints ai'e 
shown as relative probabilities in each panel. The IGM temperature inferred 
from the wavelet amplitude PDF (dot dashed, blue) is in broad agreement 
with the flux PDF value (dashed red). The cosmological parameters are un- 
constrained, as is Tgg for the wavelet amplitude PDF. The joint constraints 
are also shown (solid, black). 

In Figure[TT]we attempt an analogous comparison of our con- 
straints (thick black error bars) with the available models for 7(2) 
as well as previous measurements. Our resul ts are shown with 
the analytical reionisation models of iHui &"G nedin (1997) (with 
Treion = 25 000 K at z = 6) and the extended He ii reionisation 
m odels 'LI' and ' Lib' from the radiative transfer simulations 
of iMcOuinn et al.1 ( |2009|) , which implement 7 = 1.3 and 1.0 
at z = 6, res pectively. We compare our results with the lim- 
its from | Lidz et al. ( 2010). and measurements from [ Ricotti et al.l 
( l2000l) . ISchave et alj ( I2OOT) and lMcDonald et al.l jloOlh . Overall, 
it is clear that 7 remains relatively poorly constrained, although 
the data appear to prefer a temperature-density relation which is 
shallower than the 7 ^ 1.6 expected if Hen reionisation occurred 



• We have checked that the impact of a continuum which has 
been placed 3 per cent too low on the wavelet and flux PDF is 
small compared to the effect of varying other free parameters 
such as To, 7 and T^g. The flux PDF is indeed more sensitive to 
changes in the continuum placement, but the effect remains small 
and is minimal within the flux range of F = [0.1, 0.8] we fit in our 
analysis. We conclude it is unlikely that the continuum placement 
is fully responsible for the systematic offset found in the wavelet 
and flux PDF constraints. 

• We have explicitly confirmed that varying cosmological 
parameters within a narrow range has little impact on the wavelet 
amplitude PDF, with the strongest effect being that of ag in our 
lowest redshift bin. We also confirm that there is a significant 
degeneracy between the parameters To, 7 and Tch inferred from 
the wavelet amplitude PDF. 

• The flux PDF puts a much stronger constraint on t^s com- 
pared to the wavelet PDF. We therefore perform a joint analysis of 
the flux PDF and wavelet PDF in order to add the reft constraint 
derived from the flux PDF. We find the joint constraints on the IGM 
temperature at mean density. To, obtained at 2 = [2.1, 2.5, 2.9] 
are in good agreement with other recent measurements. The con- 
straints are consistent with the models of McOuinn et al. (2009i), 
in which an extended He 11 reionisation epoch completes around 
2 = 3, driven by quasars with an EUV index of a ~ 1.6. We 
have also performed a rudimenta ry comparison with th e recently 
proposed blazar heating models of lPuchwein et al.l ( l201ll) . and find 
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(0.36 ± 0.09) 
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0.27 ±0.02 



Table 2. Icr constraints on the IGM thermal parameters for the wavelet PDF, flux PDF and joint analysis. We find broad consistency between the methods, 
with the flux PDF favouring a slightly lower value of 7. Entries in parentheses are prior dominated (no detection) and the limits quoted are 95% confidence. 



their weak blazar heating matches our To constraints most closely. 

• We find the slope of the temperature-density relation obtained 
from the joint analysis is consistent with 7 ~ 1.1-1.3, although the 
uncertainties on this measurement remain large and remain consis- 
tent with an inverted (7 < 1) temperature-density relation within 
1-1.5(7. A more precise measurement of the temperature-density 
relation will be necessary for stringently testing competing IGM 
heating models, such as the volumetric heating rate from blazar 
heating which p roduces a strongly inverted temperature-density re- 
lation by z = 2 iChang et all201lh . 



Overall, our results are consistent with previous observations 
that indicate there may be additional sources of heating in the IGM 
at 2: < 4. This heating could be due to a number of effects such as 
an extended epoch of He 11 re i onisation which has yet to complete 
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2011 ), heating of the low density IGM by blazars dPuchwein et al 
201 ih or feed back either in the form o f galactic winds and/or AGN 
feedback ( To matore et"al]|2010l : ISoo th & Schave 2012). We note, 
however, the latter will provide only a partial explanat ion due to the 
small volume filling factor of the shock-heated gas jTheuns et al.l 
I2OO2I) . 

We believe t hat analyses of the flux PDF ( Vi el et alJboOSl : 
ICalura et al.ll2012l) . the wavelet P DF tLidz et al . 20101) and other 
techniques such as the curvature |Becker et"ail 2011 ) and absorp- 
tion line widths dSchave et al1l200oir provide generally consistent 
constraints on the IGM thermal state at 2; < 3 when applied to 
high resolution spectra. This is encouraging given the wide range 
of different methodologies used in the existing literature. However, 
significant uncertainties on measurements of the IGM thermal state 
remain. Ideally, one should also aim at reaching full consistency 
between different data sets (high and l ow reso l ution Q SO spectra) 
as well as different methods (see e.g. I Viel et aklboOSl) . In the near 
future, large scale surveys of the Lyman-Q forest at moderate spec- 
tral resolution, s uch as the SDSS B aryon Oscillation Spectroscopic 
Survey (BOSS) dSlosar et alj|201 ll) . will provide valuable new in- 
sights into the physical state of the IGM with a high degree of sta- 
tistical precision. Fully understanding the potential systematic un- 
certainties associated with these measurements are therefore vital 
for further unravelling the thermal history of the IGM following 
reionisation. 
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Figure 10. A comparison of our constraints with literature models for T(z). 
Our constraints are the thick data points: Wavelet PDF 2cr (outer, red), joint 
analysis 2(t (inner, black). ETidz et al. (2010) wavelet analysis 2(t (trian- 
gles), and the lSecker et alj )201 lb curvature analysis 2cr (squares) which as- 
sumes 7 = 1 .3. The shaded regions (upper panel) are two He II reionisation 
models from lMcOuinn et al] j2009l) ('S4b' upper, 'Dl' lower with quasar 
UV spectral indices of 0.6 and 1.6, respectively) while the lines (lower 
pane l) show the temperatur e at mean density of the blazar heating models 
from lPuchwein et alj J201 ll) . 



12 A. Garzilli, J.S. Bolton, T.-S. Kim, S. Leach and M. Viel 



1.8 




■ ■ ■ ■ 

2 3 4 5 

Redshift, z 



Figure 11. A comparison of our constraints on the slope of the 
temperature-density relation with various models for -/{z) and other ob- 
servational constraints from the literature. Our 2a constraints and limit 
from the joint analysis (thick, bl ack error bars and arrow), can be com- 
pared with the 2a limits from Mdz et al (brown arrows), the 
1(T constraints from Ricotti et al. 1200(f) (triangles, navy), the la con- 
strain ts from Schave et al. (2000) (squares, grey) and the Icr constraints 
from iMcDonald et ali (200^) (crosses, orange). The extended solid line 
shows a model with sudden He II reionisation heating the IGM to 25 000 
K at 2 = 6 (solid), while the broken lines are the ' LI' and 'Lib' models 
of He II reionisation (solid and dotted) from iMcOuinn et alj i2009l) . which 
implement 7 = 1.3 and 1.0 at 2: = 6, respectively. 
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APPENDIX A: ADDITIONAL MATERIAL 

Al Dependence of the wavelet amplitude PDF on 
cosmological parameters 

In this appendix we explicitly show the wavelet amplitude PDF 
depends almost negligibly on the cosmological parameters in our 
analysis. Figure IaTI shows the effect of varying Ho, erg, Sim, and 
Us on the wavelet PDF. 

A2 Simulation parameters and convergence tests 

In Table lAll we list the parameters of our simulations, and cat- 
alogue how they have been used in the interpolation scheme de- 
scribed in Section |4m Using simulations Rl, R2, R3 and C15 we 
have tested the convergence of the wavelet amplitude PDF with gas 
particle mass and box size. Fi sure lA2l demonstrate s the stability of 
the wavelet amplitude PDF under this change. 

A3 Continuum test 

Fi gure |A3l demonstrates the effect of lowering the continuum on the 
simulated spectra by 3 per cent on the wavelet amplitude and flux 
PDFs. A brief discussion of the results from this test is provided in 
SectionlSJ] 
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Figure Al. The wavelet amplitude PDF obtained from simulated spectra varying a single cosmological parameter about the fiducial model. Top left: Varying 
-ffo- Top right: Varying ag. Bottom left: Varying Om- Bottom right: Varying ris. The wavelet amplitude PDF has a very weak dependence on the cosmological 
parameters, but with a mild dependence on ag in our lowest redshift bin 2 = 2.1. 



Table Al. Parameters of the simulations used for describing the thermal state of the IGM: L is the comoving box-length; A'^ the number of gas and dark matter 
particles in the simulation; Mgas is the mass of each gas particle in the simulation box; and ^ are the scaling parameters used to modify the photo-heating 
rates. To and 7 are not a priori parameters of the simulations, but are determined by fitting the temperature-density relation at each redshift. The table is 
divided into five sections: our reference model (D15); simulations for interpolating over 7; simulations for interpolating over To; the validation set; resolution 
and box-size checks. 
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Figure A2. Convergence of the wavelet amplitude PDF with increasing gas particle mass at fixed box size (upper row) and increasing box size at fixed gas 
particle mass (lower row). Our conclusions from these tests are described in Section l2!2l 





Figure A3. The effect of the continuum level on the wavelet amplitude and flux PDFs at redshift z = 2.9. In the left panels we show the wavelet PDFs for 
the simulations D15 (above) and D07 (below); in the right panels we show the flux PDFs for the same simulations. The dotted lines are the results for the 
native continuum level in the simulations (indicated with regular) whereas the dashed lines are for a continuum level which has been lowered by 3 per cent. 
For comparison, we also show with a solid line the data with the error bars used in the two analyses. 



